Inverse statistics in turbulence 
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Inversion formula of multifractal energy dissipation in 3D fully developed turbulence 
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The concept of inverse statistics in turbulence has attracted much attention in the recent years. 
It is argued that the scaling exponents of the direct structure functions and the inverse structure 
functions satisfy an inversion formula. This proposition has already been verified by numerical data 
using the shell model. However, no direct evidence was reported for experimental three dimensional 
turbulence. We propose to test the inversion formula using experimental data of three dimensional 
fully developed turbulence by considering the energy dissipation rates in stead of the usual efforts on 
the structure functions. The moments of the exit distances are shown to exhibit nice multifractality. 
The inversion formula between the direct and inverse exponents is then verified. 

PACS numbers: 47.53.+n, 05.45.Df, 02.50.Fz 



>>'. 
i 

O 

>v 
43 



> 

00 

m 
o 
\o 
o 

o 

>>: 
43 

Oh- 



X 



Intermittency of fully developed isotropic turbulence is 
well captured by highly nontrivial scaling laws in struc- 
ture functions and multifractal nature of energy dissipa- 
tion rates 0. The direct (longitudinal) structure func- 
tion of order q is defined by S q (r) = (Av\\(r) q ), where 
Au|| (r) is the longitudinal velocity difference of two po- 
sitions with a separation of r. The anomalous scaling 
properties characterized by S q (r) ~ r^ ? ) with a nonlin- 
ear scalingexponent function £(q) were uncovered exper- 
imentally 



inverse scaling exponents, 
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While the direct structure functions consider the sta- 
tistical moments of the velocity increments Av measured 
over a distance r, the inverse structure functions are con- 
cerned with the exit distance r where the velocity fluc- 
tuation exceeds the threshold Aw at minimal distance 
0| . An alternative quantity is thus introduced, denoted 
the distance structure functions or inverse structure 
functions HSl, that is, T q (Av) = (r q (Av)}. Due to the 
duality between the two methodologies, one can intu- 
itively expected that there is a power-law scaling stating 
that Tp(Av) ~ Au 0(p) , where <p(p) is a nonlinear con- 
cave function. This point is verified by the synthetic 
data from the GOY shell model of turbulence exhibit- 
ing perfect scaling dependence of the inverse structure 
functions on the velocity threshold Q . Although the in- 
verse structure functions of two-dimensional turbulence 
exhibit sound multifractal nature 5], a completely differ- 
ent result was obtained for three-dimensional turbulence, 
where an experimental signal at high Reynolds number 
was analyzed and no clear power law scaling was found in 
the exit distance structure functions 4]. Instead, differ- 
ent experiments show that the inverse structure functions 
of three-dimensional turbulence exhibit clear extended 
self-similarity 

For the classical binomial measures, Roux and Jensen 
have proved an exact relation between the direct and 



C(q) = -P 

4>(p) = -q 



(i) 



which is verified by the simulated velocity fluctuations 
from the shell model. This same relation is also derived 
intuitively in an alternative way for velocity fields [Toj ] . A 
similar derivation was given for Laplace random walks as 
well [JlJ- However, this prediction (JJJ is not confirmed 
by wind-tunnel turbulence experiments (Reynolds num- 
bers Re = 400 - 1100) Q. We argue that this dilemma 
comes from the ignoring of the fact that velocity fluc- 
tuation is not a conservative measure like the binomial 
measure. In other words, Eq. Q can not be applied to 
nonconservative multifractal measures. 

Actually, Eq. Jjl is known as the inversion formula and 
has been proved mathematically for both discontinuous 
and continuous multifractal measures [12L 1 1 3| . Let fi be 
a probability measure on [0, 1] with its integral function 
M (t) = /i([0, t]). Then its inverse measure can be defined 
by ' 



= Aft 00 = 



inf{i 
1, 



M(t) > s}, 



if s < 1 
if s = 1 



(2) 



where Aft(s) is the inverse function of M(t). If /i is self- 
similar, then the relation p, = YH=i Pi^C^i (')) holds, 
where m^s are similarity maps with scale contraction ra- 
tios Ti G (0, 1) and ^Zi—xPi = 1 with pi > 0. The multi- 
fractal spectrum of measure /x is the Legendre transform 
/(a) of r, which is defined by 
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It can be shown that 0, ^| , the inverse measure ^ is 
also self-similar with ratio r\ = pi and p\ — r^, whose 
multifractal spectrum /f(af) is the Legendre transform 
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of which is defined implicitly by 

n 

E(pJ) p (^)- e - 1 ■ (4) 

i=l 

It is easy to verify that the inversion formula holds that 
f r(g) = -p f5 s 

Two equivalent testable formulae follow immediately that 
r(q) = -0-H-q) (6) 

and 

9{p) = -t-\- P ) . (7) 

Due to the conservation nature of the measure and its 
inverse in the formulation outlined above, we figure that 
it is better to test the inverse formula in turbulence by 
considering the energy dissipation. 

Very good quality high-Reynolds turbulence data have 
been collected at the SI ONERA wind tunnel by the 
Grenoble group from LEGI Q. We use the longitudinal 
velocity data obtained from this group. The size of the 
velocity time series that we analyzed is N w 1.73 x 10 7 . 

The mean velocity of the flow is approximately (v) = 
20m/s (compressive effects are thus negligible). The root- 
mean-square velocity fluctuations is v rms = 1.7m/s, lead- 
ing to a turbulence intensity equal to / = v lms /(v) — 
0.0826. This is sufficiently small to allow for the use of 
Taylor's frozen flow hypothesis. The integral scale is ap- 
proximately 4m but is difficult to estimate precisely as 
the turbulent flow is neither isotropic nor homogeneous 
at these large scales. 

The Kolmogorov microscale rj is given by [l4j 77 = 

- 2/ %2 "I 1/4 

rg77a7arre\ = 0.195mm, where v = 1.5 x 10 _5 m 2 s _1 

15((ov i oty) 

is the kinematic viscosity of air. dv/dt is evaluated by 
its discrete approximation with a time step increment 
dt = 3.5466 x 10 _5 s corresponding to the spatial resolu- 
tion e = 0.72mm divided by (v), which is used to transform 
the data from time to space applying Taylor's frozen flow 
hypothesis. 

The Taylor scale is given by [13 A = ((a j^)°V /2 = 
16.6mm. The Taylor scale is thus about 85 times the 
Kolmogorov scale. The Taylor-scale Reynolds number 
is Re\ = = 2000. This number is actually not con- 
stant along the whole data set and fluctuates by about 
20%. 

We have checked that the standard scaling laws pre- 
viously reported in the literature are recovered with this 
time series. In particular, we have verified the validity of 
the power-law scaling E(k) ~ k^ 13 with an exponent /? 
very close to | over a range more than two decades, sim- 
ilar to Fig. 5.4 of provided by Gagne and Marchand 



on a similar data set from the same experimental group. 
Similarly, we have checked carefully the determination of 
the inertial range by combining the scaling ranges of sev- 
eral velocity structure functions (see Fig. 8.6 of 1, Fig. 
8.6]). Conservatively, we are led to a well-defined inertial 
range 60 < r/77 < 2000. 

The exit distance sequence r(SE) = {rj(5E)} for a 
given energy threshold SE can be obtained as follows. 
For a velocity time series {vi = v(ti) : i = 1,2, •••}, 
the energy dissipation rate series is constructed as {Ei = 
(vi+i — Vi) 2 }. We assume that Ei is distributed uniformly 
on the interval [U, t i+ i). A right continuous energy den- 
sity function is constructed such that e(t) = Ei for t G 
[ti,t i+ i). The exit distance sequence {rj(SE)} is deter- 
mined successively by Yjb=i r k/(v) = inf{t : L e(t)dt ^ 
j ■ 5E}. Since energy is conservative, we have 

r j (2SE) = r 2j -im + r 2j (SE) . (8) 

With this relation, we can reduce the computational time 
significantly. In order to determine ri(8E), we choose a 
minimal threshold E m m, one tenth of the mean of {Ei}, 
and obtain ri(E m i n ). Then other sequences of rt for inte- 
ger 8E/E m i n can be easily determined with relation |(HJ|. 

In Fig. ^ is shown the empirical probability density 
functions (pdf's) of exit distance r/77 f° r energy incre- 
ments E mul , 2E mul , and 4i? m i n . At a fist glance, the 
probability density functions are roughly Gaussian, as 
shown by the continuous curves in Fig. ^ The value of 
fio is the fitted parameter of the mean /1 in the Gaussian 
distribution. For r/77 < MOj the three empirical pdf's col- 
lapse to a single curve. However, for large r/77 > MO) the 
three empirical pdf's differ from each other, especially in 
the right-hand-side tail distributions. This discrepancy 
is the cause of the occurrence of multifractal behavior of 
exit distance, which we shall show below. 

An intriguing feature in the empirical pdf is emergence 
of small peaks observed at r/e = 1, 2, • • ■ in the tail dis- 
tributions. Comparably, the pdf of exit distance of multi- 
nomial measure exhibits clear singular peaks. Therefore, 
these small peaks in Fig.^can be interpreted as finite-size 
truncations of singular distributions, showing the under- 
lying singularity of the dissipation energy, which is con- 
sistent with the multifractal nature of the exit distance 
of dissipation energy. 

According to the empirical probability density func- 
tions, the moments of exit distance exist for both posi- 
tive and negative orders. Figure illustrates the double 
logarithmic plots of [T p (r/ r)] 1 /^ 1 ) versus SE/E for 
different values of p. For all values of p, the power-law 
dependence is evident. The straight lines are best fits to 
the data, whose slopes are estimates of 9(p)/(p — 1). 

The inverse scaling exponent 9(p) is plotted as trian- 
gles in Fig. against order p, while the direct scaling 
exponent r(q) is shown as open circles. We can obtain 
the function — r _1 (— p) numerically from the r(g) curve, 
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which is shown as a dashed line. One can observe that 
the two functions 0(p) and — r _1 (— p) coincide remark- 
ably, which verifies the inverse formulation JSJ. Similarly, 
we obtained —9~ 1 (—q) numerically from the 9(p) curve, 
shown as a solid line. Again, a nice agreements between 
r(g) and —6~ 1 (—q) is observed, which verifies JJJ. 

In summary, we have suggested to test the inversion 
formula in three dimensional fully developed turbulence 
by considering the energy dissipation rates in stead of the 
usual efforts on the structure functions. The moments of 
the exit distances exhibit nice multifractality. We have 
verified the inversion formula between the direct and in- 
verse exponents. 
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are grateful to J. Delour and J.-F. Muzy for help in pre- 
processing these data. This work was partly supported 
by the National Basic Research Program of China (No. 
2004CB217703) and the Project Sponsored by the Sci- 
entific Research Foundation for the Returned Overseas 
Chinese Scholars, State Education Ministry. 
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FIG. 1: (Color online) Empirical probability density functions 
of exit distance r/r) for energy increments SE = iJmin, 2_E m i n , 
and 4£" m i n . The value of fio is the fitted parameter of /j, in 
the Gaussian distribution. 




FIG. 2: Double logarithmic plots of [T p (r/ £ r)] 1/(p_1) versus 
SE/E for different values of p. The straight lines are best fits 
to the data. 




p,q 



FIG. 3: (Color online) Testing the inversion formula of tur- 
bulence dissipation energy. 



